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We develop the method of calculation of the leading chiral (infrared) logarithms to an arbitrary loop order 
for various form factors of Nambu-Goldstone bosons. The method is illustrated on example of scalar and 
vector form factors in massless 4D 0{N + 1)/0{N) cr-model. The analytical properties of the form factors are 
derived. The leading chiral (infrared) logarithms are summed up in the large limit. 
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Recently we developed a new puissant method [T] 
for calculations of leading chiral (infrared) logarithms 
in a wide class of non-rcnormalizable massless field 
theories. The method has been applied to the am- 
plitude of Nambu-Goldstone boson scattering in 4D 
0{N + 1)/0{N) sigma-model defined by the following 
Lagrangian: 



(1) 



'where the fields are constrained by the relation + 
^^j^ tt'^tt" — . The amplitude of tttt scattering can 
,be decomposed into the invariant tensors of 0{N) group 
'as follows: 
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The amplitudes A, B and C are functions of the Mandel- 
stam variable s and cm. scattering angle B. In Ref. [T] 
we derived these amplitudes in the leading logarithms 
(LLs) approximation. The LLs approximation consists 
in the summation of contributions of the type ~ [s In s]" 
in the low-energy expansion of the amplitude. Such con- 
tributions arise from the n-loop Feynman graphs of the 
effective theory ([1]). At the first glance, a mission is 
impossible - to compute the n-loop graphs in a non- 
renormalizable theory. However, in the LLs approxima- 
tion this task can be accomplished as it can be reduced 
to the calculations of the spectrum of anomalous dimen- 
sions of the 0{N) symmetric composite operators made 
of four pion fields in an free field theory The result 
Ref. [T] for the tttt scattering amplitudes in LLs approx- 
imation can be presented in the form of the partial wave 
decomposition as: 
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We introduced dimensionless expansion parameter S = 
s/{4:TtF)'^ and we collected only contribution with max- 
imal power of chiral logarithm L = ln(— /i^/s), Pi{cos9) 
are Legendre polynomials. The crossing matrix fl^^ has 
the following form: 
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Pi{x)ix-ir. (4) 



The LLs coefficients tUnj satisfy non-linear recursion re- 
lation [1]: 
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which allows us to express the higher coefficients w„j 
through the coefficients with lower values of n, start- 
ing with Woo = 1- [We remind that n enumerates the 
loop order and j < {n -\- Mod(n, 2))]. The coefHcients 
B) ' ' ' are given by: 
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The recursive relation ([5]) allows a very fast computation 
of LLs. For example, the 33-loop chiral LL is computed 
in a dozen of seconds on a PC ^\ 

In present paper we develop the general method for 
calculation of LL's corrections to the form factors of the 
Nambu-Goldstone bosons (pions). We calculate here 
LL's for the scalar and vector form factors of pions in 
the massless 0{N + 1)/0{N) cr-model (P). This model 
for iV = 3 is equivalent to the chiral SU{2) x SU{2) 
model which describes the leading low-energy interac- 
tion of Nambu-Goldstone bosons (pions) of QCD in the 
chiral limit [5]. The chirally odd scalar form factor in 
the effective theory ([1]) is defined as: 



{0\J{0)\tt^{pi)7t''{p2)) = S'^' 2BFs{W^), 



(7) 



where the scalar form factor Fs{W'^) depends on the 
invariant mass of two pions — (pi + ^2)^ and the 
chirally odd scalar operator J{x) is defined as: 



N 



J{x)^2BF{F -a{x)) = B^ti''{x)-k''{x) + 0{-k'^).{^) 

a=l 

Here the constant B is proportional to the order param- 
eter of spontaneously broken symmetry; in the case of 
strong interactions it is proportional to the quark con- 
densate B = - {Q\uu\Q)lF'^. 

The chirally even vector form factor in the effective 
theory ^ is defined as: 

(0|jjr''l(0)k=(pi)7r'^(p2)) = i{P2 -Pi), (9) 

where the vector current J^/^^\x) is defined as follows: 

= ^''{x)d^^\x) - ^\x)d^^''{x). (10) 

This current is nothing but the Noether current cor- 
responding to the global 0{N) symmetry of the La- 
grangian ([1]). The low energy expansion of the scalar 
FsiW"^) and the vector FyiW^) form factors has the 
following structure: 
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where /^^ are the coefficients of the low-energy expan- 
sion in the powers of dimensionless variables 



" (4^' ^^l^-W^l ' 



^'Mathematica notebook for computing LLs is available at 
|http : / / www . tp2 .rub.de / ~maximp /research/research . html | 



where /i^ is the normalization scale. 

S V 

The lowest coefficients /qq = 1 is obtained from 
tree-level calculations of the matrix elements ([7]) and 
The calculation of higher order coefficients re- 
quires consideration of the loop diagrams in effective 
theory. Generically, the coefficients can be ob- 

tained from the calculation of diagrams with number of 
loops < n with inclusion of vertices from higher order 
effective Lagrangians which contain 2r derivatives 
with r = n + 1 — k. The calculation of the LL coeffi- 
cients /^'^ = fni^ is reduced to calculation of n-loop 
diagrams with vertices generated by the the leading ef- 
fective Lagrangian ([1]). Presently, record calculations 
of the LLs coefficients for the scalar form factor (for 

= 3) are performed by Bissegger and Fuhrer [7] to 
the four- loop order with the result: 
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/f = 1, fl = fl = TTT^, l! = = ■ (12) 
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The vector form factor (for iV = 3) is know to the two- 
loop order |^: 



(13) 



Now we present a general method, which allows us 
to perform the calculation of LLs for form factors to an 
unlimited order and for arbitrary N. We present the 
method for the scalar form factor FsiW"^), therefore we 
do not write (super)subscripts S in order to simplify 
notations. 

The UV divergencies in a n-loop diagram are re- 
moved by the subtraction of lower-loop graphs with in- 
sertion of the local counterterms corresponding to the 
subdivergencies of the original n-loop diagram. See de- 
tailed discussion of the structure of the subtractions in 
Refs. [SJITD]. The local counterterms relevant for our cal- 
culations renormalize the couplings g„j of the all-order 
Lagrangian (see Eq. (5)of Ref. 1 ). After subtraction 
of the UV divergencies the low-energy expansion of the 
form factor has the following structure: 



F{W^ 



n=0 



r„(M) + /«fc(g' ^)^'" ' (14) 
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where r„ (fi) are renormalized tree level subtraction con- 
stants that depend on the renormalization scale fi. By 
g we denote the infinite set of the constants gnjip) of 
the all-order chiral Lagrangian (Eq. (5) of Ref. [IJ), by 
r we denote the tree level subtraction constants r„(/i). 
The form factor F{W'^) should be /i independent, i.e. 
^J■^-;^F{W^) = 0. This requirement leads to the fol- 
lowing equation for the coefficients fnk- 
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7„ + (G + i/)/„o = 0, 

{G + H)fnk + ik+l)fnk+l=0, 



(15) 



where the 7- functions are defined as 7n(g, r) = 
/x^^r„(^) and G and H stand for the following 
operators: 



oc 00 n o 
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with /3-functions defined as Pnjis) = tJ'^-^9nj{lJ')- 
These /3- functions were discussed in details in Ref. [T]. 
The equation p5|) has an obvious solution: 



(17) 



with lowest constant /oo = 1 fixed by the tree order cal- 
culation of the form factor with the leading Lagrangian 
([J). We see from the solution (fT7|) that, in order to 
obtain the LLs (constants /„ = /„„), we have to apply 
the operator (G + H) n times to the 7- function 7„(g, r). 
This at first glance formidable problem can be solved if 
one notes that the operators G and H act as a contrac- 
tion mapping on the space of constants g, r. As it was 
shown in Ref. [1 the operator H possesses the follow- 
ing crucial property: /f" gmj = 0, if m < n , which 
implies that H"^gn+ij — nl cj„j with constants Unj satis- 
fying the non-linear recursion ([5]) . Analogously one can 
show that the operator G possesses the "contraction" 
property G^Vm = Oifm<n — 1, which implies that the 
application of the operator G n-times to 7„ has a fixed 
point. Thus, we can write: 



G"7n(g,r)=G' 



(18) 



where Vn are the constants that are independent of cou- 
plings g, r, which we use in order to find with help of 
Eq. p7|) the coefficients in front of LLs in form factor 
expansion /„ = /„„ = t;„. For computing the con- 
stants Vn (that is equivalent to LL approximation) due 
to the "contraction" properties of the operators G and 
H only the quadratic part of the 7n(g, r) contribute to 
the Eq. (fT8|) , so that we can represent the 7- function as 
following: 



n — 1 n— m 

7„(g,r)=^ ^r("- 



■j) 



'^7n9{n — m)j • 



(19) 
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The quadratic dependence of the relevant piece of the 
7-function means that the coefficients F'^''''^ can be com- 
puted from the one loop diagram shown in Fig. 1. The 
result of the calculation for the scalar form factor is: 




Fig.l One loop diagram contributing to the 7- functions 
coefficients (|20l2ip . Filled circle denotes the countert- 
erms for the form factor and the filled square denotes 
counterterms {pi) introduced in Eq. (5) of Ref. [l]. 
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and for the vector form factor we obtain: 
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In both equations matrices Vl^^ are given by Eq. 

Now with help of Eq. (fT8|) and Eq. p9|) we can easily 
obtain recursive equations for the LL coefficients for the 
scalar and vector form factors: 



^ n— 1 n—m 

- X! X! rs"y™'^Vm^'*^(n-m-l)i> (22) 

m=0 i=o 



S V 

with /q ' =1. The coefficients of 7- functions are 
given by Eq. ([20|) and Eq. ([2T|) for the scalar and vec- 
tor form factors correspondingly. The LL coefficients 
of the Nambu-Goldstone scattering amplitude ujpi can 
be found from the solution of the non-linear recursion 
relation ([5]). 

The equation together with ([5]) provides a very 
powerful tool for calculation of the LL coefficients in 
the expansion of the form factors. The results for the 
first 6 loops for the scalar and vector form factors in the 
0{N + 1)1 0{N) sigma-model ^ are presented in Ta- 
ble 1 and Table 2 correspondingly. For iV = 3 the results 
coincide with the laborious four-loop calculation of the 
scalar form-factors obtained in Ref. [7] (see, Eq. (fT2|)) 
and with 2-loop calculations of Ref. [6] (see Eq. (fT3l) . 

In the large N limit the 0{N+l)/0{N) sigma model 
can be solved by the semiclassical methods [TT]. In 
Ref. [1] we showed that in this limit the non-linear re- 
cursion equation ([5]) has the following solution : 



5,0 1 + 



(23) 



Substituting this solution into Eq. ((22|) we can easily 
solve corresponding recursion relation in the large N 
limit with the result in > 1): 
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With these results we can perform the summation of the 
LL contributions in the large N limit. The form factors 
in these approximation have the following form: 
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■ e 
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(25) 
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Here we introduced a following short hand notation 
e = 2^7r^p 111 (~ w^) ■ scalar form fac- 

tor in the large N limit and LL approximation possesses 
a pole at e = 1 (which is actually outside the applica- 
bility of the LL approximation). This pole corresponds 
to the contribution of the auxiliary scalar field, which is 
introduced in order to solve the model ([1]) at ^ cxd 
(see e.g. [H]). The vector form factor also possesses 
a weaker singularity at e = 1, that corresponds to the 
threshold of 27r-|-auxiliary scalar field production. We 
see that the form factors in the LL approximation and 
in the limit of large number of Nambu-Goldstone boson 
A'^ have correct analytical properties. Let us demon- 
strate that in general case the form factors obtained in 
LL approximation with help of our recursion equation 
satisfy all constraints imposed by analyticity on the 
form factors. 

Now we present general solution of recursion rela- 
tion (|22p that expresses LL coefficients for form factors 
f^'^ in terms of the LL coefficients of tttt scattering 
amplitude w„.j. To this end we introduce two types of 
generating functions for the coefficients: 



(26) 
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which satisfy obvious conditions Ws.v{0) = 1 and 
^s,v{0) — 1. The scalar and vector form factors in the 
LL approximation can be express in terms of generating 
functions v as follows: 



(27) 



The recursion relation ([22|) can be reduced to the dif- 
ferential equation for the generating functions ([26]) that 
has the following solution: 



^s.y{x) = exp 



dy Wsyiy) 



(28) 



From Eq. ^ we can conclude that the lowest / = 0, f 
partial wave amplitudes tj (s) — ^ ^e^**!^('^) — 1^ can be 



expressed in terms of the generating functions Wsy{x) 
in the following way: 

tO(s) = nSWsiSL), tl{s) = TrSWv{SL) . (29) 

[For notations of kinematical variables see text after 
Eq. (O.] We see that due to the relations (|27l28p the 
form factors in the LL approximation can be expressed 
in terms of partial waves for the tttt scattering in the 
following way: 




In 
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Noting the following small asymptotic to the lead- 
ing logarithms accuracy: 



/•OO 

lim / ds 
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we can put representation of the form factors in terms of 
the TTTT partial wave amplitudes (|30p in the form which 
corresponds to the Omnes solution [12] of the dispersion 
relations for the form factors: 



M/2 



(31) 



This demonstrates that the solution of the recursion re- 
lation provides the form factors with correct ana- 
lytical properties. 

In summary, we suggested new method to compute 
leading infrared logarithms to arbitrary loop order for 
the form factors of the Nambu-Goldstone bosons of an 
effective field theory in four dimensions. The method is 
demonstrated on the example of the scalar and the vec- 
tor form factors in 4D 0{N + 1)/0{N) cr-model. The 
proposed method can be applied to a wide range of ef- 
fective theories. 
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Table 1: LL coefficients for the scalar form factor 
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fiiN = 3) 







1 
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N 1 
2 2 
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43 
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N'-' 29 N 1 11 
4 72 ^ 72 
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^3 271Af^ 1 7N 89 
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Af" 6382513Af^ i 3785803199Af* 7206506437Af^ i 11173397867Ar^ 630301337Af i 255705409 


428652000 


64 84672000 ' 22861440000 34292160000 ' 68584320000 8573040000 ' 17146080000 


Table 2: LL coefficients for the vector form factor 
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1 

6 


1 

6 




2 
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N 1 
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240 1200 ' 1555200 311040 
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Af^ 20917Af'' 1 10524683A''^ 20818337V , 2729 




163296000 


672 4536000 ' 1632960000 544320000 ' 2551500 
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937784623 


N" 24313A'* 1 8249471W 478853567Ar^ , 34922611Ar 1164883001 




41150592000 


1792 9408000 ' 1524096000 68584320000 ' 6531840000 411505920000 





1. N. Kivel, M. V. Polyakov and A. Vladimirov, Phys. Rev. 
Lett. 101 (2008) 262001 ;arXiv:0809.3236 [hep-ph]]. 

2. S. Weinberg, Physica A 96 (1979) 327. 

3. S. Weinberg, Phys. Rev. 166 (1968) 1568. 

4. P. Langacker, H. Pagels, Phys. Rev. D 8 (1973) 4595. 

5. J. Gasser, H. Leutwyler, Annals Phys. 158 (1984) 142. 

6. J. Bijnens, G. Colangelo, G. Ecker, J. Gasser and 
M. E. Sainio, Nucl. Phys. B 508 (1997) 263 [Erratum- 
ibid. B 517 (1998) 639] arXiv:hep-ph/970729i; . 

7. M. Bissegger and A. Fuhrer, Phys. Lett. B 646, 72 
(2007) [arXiv:hep-ph/0612096| . 

8. J. Bijnens, G. Colangelo and P. Talavera, JHEP 9805 
(1998) 014 arXiv:hep-ph/980538 9 . 

9. M. Buchler and G. Colangelo, Eur. Phys. J. C 32, 427 
(2003) arXiv:hep-ph70309049 . 

10. N. Kivel and M. V. Polyakov, Phys. Lett. B 664 (2008) 
64 [arXiv:0707.2208 ll [hep-ph]]. 

N. Kivel, M. V. Polyakov and A. Vladimirov, Phys. Rev. 
D 79 (2009) 014028 ;arXiv:0809.2064 [hep-ph]]. 

11. S. R. Coleman, R. Jackiw and H. D. Politzer, Phys. Rev. 
D 10 (1974) 2491. 

12. R. Omnes, Nuovo Cim. 8 (1958) 316. 



